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Abstract 

This paper presents estimations on the possibility of detection of ellipsoidal variations by means of 
measuring brightness of the star distorted by a close massive planet using Wilson-Devinney method. 
The problem was already discussed by Phafl et al. (2008) and earlier by Loeb & Gaudi (2003). The 
effect is well known in the case of binary stars where it can produce light curves with amplitutudes 
of ellipsoidal variations of about 0.1 mag for distorted stars. For planets the effect is very small and 
usually less than 0.0001 mag. The detection of an exoplanet, by searching for small amplitude ellipsoidal 
variations, will be very difficult and affected by other photometric effects; however, it maybe possible 
for some extreme cases. Observations of ellipsoidal variations can provide additional constraints on the 
model of the system. Light curves for few star/planet systems have been calculated using PHOEBE 
eclipsing binary software based on Wilson-Devinney method. As an example of ellipsoidal variations the 
synthetic light curve for rBootis is presented. The amplitude of ellipsoidal variation is 0.01 mmag. The 
companion is massive (7.3 Mj up ) and short-period hot Jupiter. 

1 Introduction 

During the last decade we have observed a very fast development of extrasolar planet research. The new 
methods of detection have been successfully used. The observational techniques give us the possibility 
to detect more and more shallow effects like transits of a planet in front of the hosting star, a direct 
imaging of planets and measuring the chemical abundance in the planet atmosphere. In this paper I 
will discuss the question whether we can observe light curve variations of the hosting star distorted by 
the gravity field of a massive planet for some of the already known systems. The detailed analytical 
formulation of ellipsoidicity effect caused by substellar companion is given by Pfahl et al. (2008). A 
close massive companion of the star can produce tides analogue to the ones present in close binary stars. 
The effect is of course much weaker and difficult for detection. Actually, the presence of ellipsoidal 
variations among the transit candidates is treated as an evidence of a non planetary behaviour of the 
eclipses. Patzold et. al (2004) describe orbital period changes in planetary systems caused by gravity 
tides induced by exoplanets. Some short-period exoplanets, for example rBootis, are probably in tidal 
lock that demonstrates the significance of the tidal interaction. We can model ellipsoidal variations using 
a code based on Wilson-Devinney (WD) method developed for eclipsing binary stars. This method has 
already been used for modeling eclipses and radial velocity curves of transiting planets - one of them 
being OGLE-TR-56 (Voccaro and Van Hamme 2005). 

2 Roche lobe and the massive planet 

I will use a standard Wilson-Devinney method (PHOEBE software, Prsa and Zwitter 2005) to describe 
the distortion of the hosting star by its planetary companion. In the case of eclipsing binaries we have a 
typical mass ratio of about q = 7712/mi ~ 1. For planetary systems, where mi >> m,2 we usually have a 
mass ratio less than 0.01. The ellipsoidal light curve variations in binary systems are dependent on the 
level of asymmetry of the distorted component. We have higher amplitude of light variations for stars 
that significantly fill their Roche lobes. For mass ratio of about 1 we observe distorted Roche lobes with 
high asymmetry near Lagrangean point L%. For q < 0.01 the lobe of the star is almost spherical (for slow 
rotation) with a small asymmetry near L\. The shape of the Roche lobe also depends on a rotational rate 
of the components. In WD method we use the synchronicity parameter /, defined as the ratio between 
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Figure 1: Synthetic light curve of the hypothetical system for two eccentricities, inclination of 79° and 
synchronous rotation. The gray dots present the numerical noise for the circular orbit. 



the axial and orbital revolution period. For high values of this parameter, / ~ 1 (synchronous rotation), 
we have a flattened star - this implies higher asymmetry because of smaller distance between the stellar 
surface to L\. In table[T]we have listed values of amplitudes of ellipsoidal variations for two synchronicity 
parameters / of both the hosting star that fills its Roche lobe and for the MS star. The amplitude also 
depends on the gravity and limb darkening coefficients. The main contribution in ellipsoidal variation 
amplitude is connected with the stellar shape (> 90%) and gravity brightening (< 10%). We can expect 
notable ellipsoidal variations in two cases, considerable mass ratio or high filling of the Roche lobe. The 
distortion computed by the Roche model is slightly overestimated for evolved stars because the mass is 
not concentrated in the barycenter but spread in a distorted star. 

Let us take a system with a relatively distorted star characterized by a mass ratio of 0.001, a short 
orbital period of Id and a semi-major axis of 5 R©. The main parameters of the system are close to 
those of OGLE-TR-56. The mass of the planet is 0.0013 M (1.3 Mj up ). The host star is slightly evolved 
with a mass of 1 Mq and a radii of 1.3 Rq. The obtained light curve has peak to peak amplitude of 
about 0.00006 mag. The distorted star has the following dimensions: i? s i<j c = 1.32179, -Rpoint = 1.32184, 
Rback = 1.32183, -Rpoio = 1.30846 in solar units. Two curves are presented, one for circular and one for 
eccentric orbit (fig. [1]). For a more evolved star with radii 2.5 R© we have peak to peak amplitude of 
0.0005mag. 

We have two minima of different depths, which is the result of different surface temperature of the 
star near Lagrangean points Li (lower) and L2 (higher). Two maxima are different because of the 
eccentric orbit that gives a variable separation of components for different phases. This implies changes 
of radii/temperature of the star. We can also see another effect of the orbital eccentricity - shifted minima 
and maxima. The reflection effect was switched off, just as in the other modeled objects. The numerical 
noise in LC is on the level of 10 -5 for stellar surface grid raster n B t ar = 60 (maximal for PHOEBE). For 
the future aplications the exoplanet mode with denser grid can be of use for modeling shallow effects. 

3 Ellipsoidal variations of known systems 

Among the discovered systems listed in exoplanet. eu data-base some stars have been checked for ellip- 
soidal variations using WD simulations. The star with higher amplitude is HD 41004B. For the generation 
of its light curve the results of Zucker et al. (2004) were used. Around the host star there is a massive 
companion (18.5 Mj up ) orbiting, the mass of which is above the limit of 13 Mj up usually applied for 
planets; therefore, we can classify this object as a brown dwarf. The semi major axis of the orbit is very 
small 0.018 AU and the orbital period is 1.3d. Obtained LC presents variations with peak to peak am- 
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Figure 2: Both components of the rBootis system inside their Roche lobes for inclination 80°, which presents 
better the flatness of the lobes. 



Table 1: Maximal V amplitude of ellipsoidal variations vs. mass ratio for the star that fill its Roche lobe 
and for the main sequence star. Inclination is 90°, component 2 is small and spherical - does not involve the 
curve shape. The host star has one solar mass, the semi-major axis is 5 R© and the period ~ 1.3d. 



mass ratio 


evolved 


main sequence 


q 


f=0.1 


f=l 


f=0.1 


f=l 


1 


0.27 


0.28 


0.022 


0.022 


0.1 


0.12 


0.13 


0.0025 


0.0025 


0.01 


0.025 


0.044 


0.00025 


0.00025 


0.001 


0.0045 


0.013 


10~ 5 


10~ 5 


0.0001 


0.0004 


0.0032 







plitude of 0.0002 mag (V). The theoretical amplitudes for other examined stars are: rBootis - 0.00001, 
OGLE-TR-56 - 0.00006 mag. 

The former - rBootis - is a well studied system discovered in 1997 (Butler et al. 1997). For the 
WD simulations system parameters (table [2} listed in Leigh et al. (2003) were adopted. This work 
presents spectroscopic search for reflected light. The projected orbital velocity amplitude was found to 
be K p — 97 ± 8fcms _1 (K s = 469 ± 5ms" 1 ) The most probable orbital inclination is i = 37 ± 5°, which 
implies planetary mass M p — 7.28T0.83 Mj up . The synthetic light curve for two inclinations is presented 
in figure[3] The higher possible peak to peak amplitude of ellipsoidal variations without eclipses (i = 80°) 
is 3 ■ 10 _5 mag and for inclination 37° we obtain 1 • 10~ mag variability. The eccentricity is small and 
unnoticeable on the light curve. The shape of the star depends on the filling of the Roche lobe (fig. [2). 
During the evolution rBootis will expand and fill equipotential surfaces with higher asymmetry, which 
will increase the ellipsoidal variations (table 

Future mass photometry space missions, for example Kepler (Christensen-Dalsgaard et al. 2006), will 
be able to detect such level of variability. The prospects of detection to be performed by Kepler mission 
have been discussed by Pfahl et. al. (2008). The differential photometry will reach precision of 20ppm for 
12mag G2V star. The detection of ellipsoidal variations will be affected by the other photometric effects 
induced by planets - chromospheric activity, scattered light and, not connected with the companion, for 
example spots and convection (the intrinsic stellar noise). The MOST satellite has detected a variability 
on r Bootis A probably induced by the planet (Walker et al. 2008). The active region close to the 
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Figure 3: Synthetic light curve of rBootis A for two orbital inclinations. The curves have been shifted for 
better comparison and smoothed to remove the " numerical noise" . 

sub-planetary point has been found. This variable spot produces brightness variations with peak to peak 
amplitude of 1 mmag. This effect is higher than ellipsoidal variations so their detection will probably 
be impossible if the active region is stable and present in the future. A similar effect was found on HD 
179949 by Shkolnik et al. (2005). If such spots are typical for 51 Peg systems, searching for ellipsoidal 
variations of hosting stars will be more difficult. The predicted level of scattered light for rBootis by 
Green et al. (2003) is two orders lower 1 • 10 mag and comparable with the ellipsoidal variations effect. 

4 Conclusions 

The future mass photometry missions can be used for detection of exoplanet candidates by searching for 
effects weaker than transits like induced spots, ellipsoidal variations or scattered light. The WD code is 
already used for modeling transits of exoplanets. This method can be usful for calculating weaker effects 
like ellipsoidal variations or spots on the host star induced by the planetary companion, and maybe 
applied for the scattered light after some modification of the code. The method can also be useful for 
modeling the Rossiter-McLaughlin effect in radial velocity curves. 
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Table 2: Parameters of rBootis system used for the synthetic curve and obtained dimensions of the compo- 
nents for double tidal lock. 



Parameter 


star 


planet 


i 




37° 


1 




0.005 




10.51 (0.0489 AU) 


Period (d) 




Q Q 1 O/l C 

o.olz4o 


e 




u.uio 






65° 


o 


7.12 


1.93 




0.54 




Mass(M ) 


1.42 


0.007 (7.28 M Jup ) 


Radii(R ) 


1.48 


0.12 (1.2 Rj up ) 


T cff (K) 


6360 


1576 (estimated) 


-Rside(R©) 


1.480000 


0.122414 


-Rpoint(R©) 


1.480034 


0.122472 


-Rback(R©) 


1.480028 


0.122472 


i?pole(R©) 


1.477926 


0.122393 



Table 3: Relation between the filling of the Roche lobe by the host star (fi) and the ellipsoidal variations 
(Am) for t Bootis. The dimensions of the star (point & side) are given in units of semi-major axis. 



Q 


radii R Q 


Am 


-^point 


-^sidc 


7 


1.5 


1 ■ 10~ 5 


0.143 


0.143 


3 


3.6 


2.5- 10~ 4 


0.341 


0.340 


2 


5.7 


1.2- 10~ 3 


0.542 


0.541 


1.6 


8.0 


0.005 


0.789 


0.771 
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